ABSTRACT
Objective:To explore the immune infiltration cells in rheumatoid arthritis (RA) synovial lesions, and to provide new research directions and therapeutic targets for the pathogenesis and treatment of RA.Methods:The three gene expression data sets GSE77298, GSE55457 and GSE1919 were downloaded from gene expression omnibus (GEO) (http://www.ncbi.nlm.nih.gov/geo), and the data were merged with Perl. The "limma" package was used to adjust batch differences. In R, "CIBERSORT" software was used to obtain the expression matrix of 22 kinds of immune cells corresponding to RA synovial tissue samples and normal synovial tissue samples were analyzed with the three packages of "e1071", "parallel" and "preprocessCore". Perl was used to screen samples with P<0.05 in the immune cell matrix. R's "barplot" function was analyzed by the percentage of 22 immune cells in samples with P<0.05. The "pheatmap" package of R was used to visualize heatmaps, and "corrplot" package was used to draw correlation heatmaps. The "vioplot" package of R was used to draw violin plots of differences via the wilcox test. Results:The results of immune cell infiltration analysis showed that in RA synovial tissue samples and normal synovial tissue samples at P<0.05, B cells naive and natural killer cells resting were under-expressed in RA synovial tissue, and plasma cells, mast cells resting, macrophages M1, B cells memory and T cells regulatory were highly expressed in RA synovial tissue. This study also found that in the same sample, the correlation coefficient between natural killer cells resting and neutrophils ( r=0.91) was the highest, indicating synergistic effect between the two. In the same sample, the correlation coefficient between macrophages M0 and plasma cells ( r=-0.88) was the lowest, indicating antagonistic effect between the two. Conclusion:The immune infiltrating cells in RA synovial lesions discovered in this study provide a certain theoretical basis and research direction for the research on the disease mechanism and treatment of RA.
ABSTRACT
Objective:To explore the potential Hub genes, key miRNAs, biological processes and related signaling pathways in the pathogenesis of osteoarthritis (OA), and provide bioinformatics basis for the pathogenesis and treatment of OA.Methods:The expression profiling chip of OA synovial tissue sample from Gene Expression Omnibus (GEO) were downloaded, differentially expressed genes (DEGs) were identified, and functional enrichment analysis was performed. A protein-protein interaction network (PPI) was constructed. STRING and Cytoscape was used for module analysis, and the Hub gene was further identified, and further miRNAs mining of the Hub gene was carried out.Results:Finally, 9 Hub genes (SOCS3, BTRC, FBXO32, KLHL22, UBE3A, HUWE1, UBR4, ANAPC5, TRIM50) and 2 key miRNAs (hsa-miR-103a-3P, hsa-miR-107) related to the progression of OA were identified .They might be potential biomarkers for the pathogenesis of OA. We also found that signal transduction, the transcriptional positive regulation of RNA polymerase Ⅱ promoter, and protein serine/threoninase activity had a certain correlation with the pathogenesis of OA. In addition, our analysis results showed that cAMP signaling pathway and Rap1 signaling pathway were also involved in the progression of OA.Conclusion:The potential biological molecules, biological processes and related pathways identified in this study may guide us for the further research on the etiology and treatment of OA.